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Abstract 

Background: Pyrethroid resistance has been slower to emerge in Anopheles arabiensis than in An. gambiae s.s and 
An. funestus and, consequently, studies are only just beginning to unravel the genes involved. Permethrin resistance 
in An. arabiensis in Lower Moshi, Tanzania has been linked to elevated levels of both P450 monooxygenases and 
(3-esterases. We have conducted a gene expression study to identify specific genes linked with metabolic resistance 
in the Lower Moshi An. arabiensis population. 

Methods: Microarray experiments employing an An. gambiae whole genome expression chip were performed on 
An. arabiensis, using interwoven loop designs. Permethrin-exposed survivors were compared to three separate 
unexposed mosquitoes from the same or a nearby population. A subsection of detoxification genes were chosen 
for subsequent quantitative real-time PCR (qRT-PCR). 

Results: Microarray analysis revealed significant over expression of 87 probes and under expression of 85 probes 
(in pairwise comparisons between permethrin survivors and unexposed sympatric and allopatric samples from Dar 
es Salaam (controls). For qRT-PCR we targeted over expressed ABC transporter genes (ABC '2060'), a glutathione-S- 
transferase, P450s and esterases. Design of efficient, specific primers was successful for ABC '2060'and two P450s 
{CYP6P3, CYP6M2). For the CYP4C16 gene, we used the primers that were previously used in a microarray study of 
An. arabiensis from Zanzibar islands. Over expression of CYP4C16 and ABC '2060' was detected though with contrasting 
patterns in pairwise comparisons between survivors and controls. CYP4G16 was only up regulated in survivors, whereas 
ABC '2060' was similar in survivors and controls but over expressed in Lower Moshi samples compared to the Dar es 
Salaam samples. Increased transcription of CYP4G16 and ABC '2060' are linked directly and indirectly respectively, with 
permethrin resistance in Lower Moshi An. arabiensis. 

Conclusions: Increased transcription of a P450 {CYP4G16) and an ABC transporter {ABC 2060) are linked directly and 
indirectly respectively, with permethrin resistance in Lower Moshi An. arabiensis. Our study provides replication of 
CYP4G16 as a candidate gene for pyrethroid resistance in An. arabiensis, although its role may not be in detoxification, 
and requires further investigation. 
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Background 

Vector control programmes employing ITNs and IRS 
are a mainstay of malaria control [1,2]. Continuous 
usage of insecticides under such malaria control pro- 
grammes and/or agricultural application of insecticides 
has resulted in the development of resistance in major 
malaria vector species [3-7]. Anopheline mosquito popu- 
lations in many parts of Africa have developed resistance 
to pyrethroids [8], which are widely used for indoor re- 
sidual spraying (IRS) and the only insecticide class avail- 
able for insecticide treated nets (ITNs). Pyrethroids have 
been shown to pose very low health risks to humans and 
other mammals, but are toxic to insects and knock them 
down (kill them), even at very low doses [9]. 

Insecticide resistance is considered a major threat to 
the continued success of IRS and ITNs, with reduction 
in efficacy of pyrethroids being a particular concern 
[8,10]. The various mechanisms by which mosquitoes 
are thought to develop resistance to insecticides include 
metabolic resistance, target-site resistance, reduced pene- 
tration and behavioural avoidance. However, metabolic re- 
sistance and target-site resistance are major mechanisms 
assumed to be responsible for insecticide resistance [11] 
and are known to contribute to pyrethroid resistance in 
malaria vectors. While target-site resistance occurs when 
the site of action of insecticide is altered such that the in- 
secticide no longer binds effectively, metabolic resistance 
involves over-expression of enzymes capable of detoxify- 
ing insecticides or amino acid substitutions within these 
enzymes, which alter the affinity of the enzyme for the in- 
secticide [12]. The most common target site resistance 
mechanism constitutes two point mutations at amino acid 
position 1014 of the voltage-gated sodium channel gene, 
resulting in either a leucine-phenylalanine (L1014F) [13], 
or a leucine-serine (L1014S) substitution [14]. 

Detoxification enzymes include glutathione-S -transferases 
(GSTs), esterases and P450 monooxygenases. GST-based 
DDT resistance by breaking down DDT to non-toxic 
products is common in a number of anopheline species 
including An. gambiae [15-17]. It has been suggested 
that GSTs may play a role in pyrethroid resistance by 
detoxifying lipid peroxidation products induced by py- 
rethroids and/or by protecting from insecticide expos- 
ure induced oxidative stress [18]. GSTs might also 
confer a secondary role in pyrethroid resistance by se- 
questering the insecticide, hence reducing the total 
in vivo concentration of insecticide [19]. More recently, 
the study of Riveron et al., 2014 [20] revealed a more 
direct role played by GSTs in conferring insecticide re- 
sistance in malaria vectors. 

Elevated levels or activity of esterase enzymes which 
hydrolyze ester bonds or sequester insecticides is one 
of the most common metabolic resistance mechanisms 
in organophosphate (OPs) resistant mosquito species. 



Malathion (OP) resistance in Anopheles culicifacies and 
Anopheles stephensi has been associated with an altered 
form of esterase that specifically metabolizes the mol- 
ecule at a much faster rate than that in susceptible 
strains [21,22]. 

Cytochrome P450-dependent monooxygenases are an 
important and diverse family of enzymes involved in the 
metabolism of numerous endogenous and exogenous com- 
pounds. Cytochrome P450 belongs to six families and 
increased transcription of genes belonging to the CYP4, 
CYP6, and CYP9 has been observed in various insecticide- 
resistant species [23]. Microarray-based approaches have 
identified three candidate P450 genes, CYP6M2, CYP6P3 
and CYP6Z2 that were found to be repeatedly over- 
produced in pyrethroid resistant populations of An. 
gambiae [24-26]. All of these genes encode for enzymes 
that are able to bind to type I and type II pyrethroids but 
only CYP6P3 and CYP6M2 were shown to metabolize the 
insecticides [27,28]. More recently, some studies demon- 
strated that CYP6M2 is also capable of metabolizing the 
organochlorine insecticide DDT in An. gambiae, hence 
demonstrating the first evidence for a metabolic cross- 
resistance in malaria vectors [29]. 

Both metabolic resistance and target-site resistance have 
been documented in An. arabiensis, one of the dominant 
vector species of malaria in sub-Saharan Africa [30,31] 
and in some areas such as the Great Rift Valley in East 
Africa, An. arabiensis is the predominant malaria vec- 
tor species [32], To date, the kdr mutations have been 
reported in several countries including Uganda [33], 
Sudan [34,35], Cameroon [36], and Tanzania [37-39]. In 
addition to target site resistance, metabolic resistance has 
also been documented in different countries [40-42]. 

The study of dynamics of insecticide resistance in Lower 
Moshi An. arabiensis from 2009 to 2013 [39] has shown 
that the population has developed resistance to all pyre- 
throids tested (permethrin, deltamethrin and lambdacyha- 
lothrin) with the presence of L1014F mutation at very 
low frequency. Tanzania has been scaling-up the use of 
pyrethroid-based LLINs which reached its universal 
coverage under Universal Coverage Campaign [43]. 
Therefore, higher pyrethroid resistance is expected fol- 
lowing higher selection pressure on this Lower Moshi 
An. arabiensis population and it is likely that pyrethroid 
resistance will spread to other areas in the country. The 
previous study by Matowo et al, 2010 [44] had revealed 
significantly higher levels of oxidase and (3-esterase en- 
zymes in wild Lower Moshi An. arabiensis than in a labora- 
tory reference strain. However, the genes responsible for 
elevated levels of these detoxification enzymes remained 
unknown, although over-expression of a few detoxification 
genes has been documented in An. arabiensis populations 
in Hai district, adjacent to Lower Moshi. [45]. Therefore, 
the aim of this study was to explore and identify specific 
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genes involved in metabolic resistance to permethrin in 
Lower Moshi populations of An. arabiensis from Lower 
MoshL 

Methods 

Mosquito sample collection and identification 

Collections were conducted in the rainy season between 
May and June 2011 in Mabogini village (37°21'E, 3°24'S), 
Mabogini is an agricultural area of about 15 km south of 
Moshi town, in which rice is cultivated in two growing sea- 
sons, the main season in which rice is sown is mid-June 
with irrigation from June to October, and a second season 
of sporadic cultivation from September to February. The 
rice paddies thus provide breeding sites for Anopheles 
arabiensis, the predominant malaria vector in the area, 
throughout much of the year. 

The larvae were collected from two specific localities, 
Mbugani (37°22' E, 3°25'S) and Harusini (37°21'E, 3°23'S) 
in Mabogini village which are 4 km. apart. Additional lar- 
val collections were made from Mbugani in August 2012. 

Anopheline larvae were also collected from temporary 
breeding habitats in Ilala (Gerezani) (39°16' E, 6°49'S) 
and Kinondoni (Kawe) (39° 13' E, 6°43'S) in Dar es sa- 
laam about 700 km from Moshi, in August and October 
2011. 

The field-collected anopheline larvae were transferred 
into an insectary and reared to adults. Emergent adults 
were split by sex before mating could occur and fed with 
10% sugar solution. 

Virgin females were preferred following the previous 
whole-body microarray experiments [46] that detected 
large transcriptional changes in mated females where the 
number of genes differentially expressed increased with 
time post mating. 

At 3 days-old, adult females from Mbugani were ex- 
posed to 0.75% permethrin for one hour using a stand- 
ard WHO tube bioassay procedure [47] with mortality 
scored 24 hours after the exposure period. Adult mos- 
quitoes from Harusini, Moshi and Dar es Salaam were 
not exposed to insecticide. At their fourth day post 
emergence, all unexposed (control) females and the sur- 
vivors (resistant mosquitoes) were killed in 75% ethanol. 
A hind leg was removed from each mosquito and stored 
over silica gel in an individually-labelled tube before 
placing the remainder of the specimens in RNA later in 
an Eppendorf tube. DNA was extracted from the legs 
that were removed from each specimen for species 
identification. 

Species identification focused on the most prevalent 
local vectors of An. gambiae sensu lato; An. arabiensis 
and An. gambiae s.s. Mosquitoes were identified to spe- 
cies level using standard PCR method [48]. Only those 
identified as An. arabiensis were further analysed in this 
study. 



Whole genome microarrays 

Experimental design 

The study was designed to comprise of two microarray 
experiments. The first involved three groups of samples 
collected in 2011 from Mbugani (insecticide exposed, 
selecting for the 24% and 54% of the most resistant fe- 
males [39] and unexposed) and Harusini (unexposed). 
Although initially suspected that Harusini might display 
a lower resistance level than Mbugani, its close proxim- 
ity made predictions ambiguous and only results from 
hybridizations of Mbugani permethrin-selected against 
Mbugani unexposed and Mbugani permethrin-selected 
against Harusini unexposed were used. In the second ex- 
periment, to investigate temporal repeatability of results 
we repeated the experiment from Mbugani (permethrin- 
selected against, unexposed) using a collection from 
2012. In this second experiment a pyrethroid-susceptible 
sample from Kinondoni and Ilala in Dar-es-Salaam [40] 
was also included, primarily to permit an allopatric com- 
parison with the population from Mbugani (insecticide 
exposed that selected for the 54% of the most resistant 
females [39]. However, we found that comparisons in- 
volving the Dar-es-Salaam sample yielded very little add- 
itional resolution for the analysis, which coupled with 
suspected technical errors in some of the hybridizations 
led to their exclusion. 

2011: Mbugani permethrin selected vs Mbugani un- 
exposed 

2011: Mbugani permethrin-selected vs. Harusini un- 
exposed 

2012: Mbugani permethrin-selected vs. Mbugani un- 
exposed 

Comparisons in experiment 1 consisted of three inde- 
pendent biological replicates while those of experiment 
2 consisted of four independent biological replicates 
(Figure 1); and two technical repeats of dye swaps to 
control for dye bias. Hybridization and analysis followed 
an interwoven loop design. 

To be judged significant a probe was required to (a) 
exhibit a false discovery rate-corrected probability of 
q < 0.05 in all three comparisons, and (b) exhibit a 
consistent directionality in expression, i.e. either Mbugani 
permethrin-selected groups > unexposed groups or 
Mbugani permethrin-selected groups < unexposed groups 
in all three comparisons. Owing to the use of within 
population (i) selected vs. (ii) unexposed comparisons, 
the magnitude of differential expression is expected to 
be relatively low, because (i) is a subset of (ii) [20] 
and therefore no fold-change threshold was applied as 
a significance criterion. 

RNA extraction and labelling 

Both four day old unfed unmated female mosquitoes 
that survived one hour exposure to 0.75% permethrin, 
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Figure 1 Interwoven microarray experimental loop design for 
experiment 1 and 2. 

V J 



held for 24 hours post-exposure and unexposed mosqui- 
toes were killed using 70% absolute ethanol and stored in 
RNAlater® at -20°C Total RNA was extracted from 
batches of ten mosquitoes using the RNAqueous^-4PCR 
Kit (Ambion) according to manufacturers instructions 
and treated with DNase I (Qiagen). The quality and quan- 
tity of the RNA was checked using a using a Nanodrop 
spectrophotometer (Nanodrop Technologies, Wilmington, 
DE, USA), and Bioanalyser (Agilent Technologies, 
Stratagene, USA) respectively. Extracted mRNA was 
dye-labelled using the Low Input Quick Amp Labelling 
Kit (Agilent Technologies, Stratagene, USA) according to 
manufacturers instructions. Labelling was repeated for 
samples that had a low cDNA yield (<15 ng/|il) and poor 
dye incorporation (<6.0 pmol/|il of each dye). 

The labelled samples were hybridized using a Gene 
Expression Hybridization Kit (Agilent Technologies). 
Washing, scanning and feature extraction were performed 
according to the manufacturer s recommendations. 

Microarray analysis 

Analysis of the microarray data was performed on cor- 
rected spot intensities using the LIMMA 2.4.1 soft- 
ware package [49] for R 2.3.0. Normalization of signal 
intensities was followed by ANOVA F-test using the 



MAANOVA package in R [50]. Microarray data are 
available in the ArrayExpress database (www.ebi.ac.uk/ 
arrayexpress) under accession number E-MTAB-2513 
(SURV vs UNEXPOSED loop design) and E-MTAB- 
2514 (SURV vs DAR loop design). 

Quantitative real time PCR (qPCR) 

Technical repeatability of microarray results for a subset 
of candidate genes was assessed using Quantitative PCR 
using the cDNA synthesised from the same RNA sam- 
ples used in the microarray experiments. The samples 
were quantified using a Nanodrop spectrophotometer 
(Nanodrop Technologies, Wilmington, DE, USA). cDNA 
was synthesised from -0.5-1 \xg of RNA using oligo(dT) 
20 (50 \iM) and Superscript III (200U) (Invitrogen) and 
purified through a DNA-binding column (Qiagen). The 
quality and quantity of cDNA was measured using the 
Nanodrop spectrophotometer (Nanodrop Technologies, 
Wilmington, DE, USA) and the cDNA samples were 
stored at -70°C until further use. 

Primers for qPCR were designed using NCBI primer 
BLAST (http://www.ncbi.nlm.nih.gov/tools/primer-blast/) 
by using Xm codes from Vector Base. Serial dilutions 
of cDNA were used for standard curve production to 
determine PCR efficiency and specificity of the primer 
pairs. The qPCR reactions were performed using the 
Agilent MXPro Real-Time PCR detection system (Agilent 
Technologies, Stratagene, USA). 

A total volume of 20 \A contained 10 \A Brilliant III 
SYBR Green, 300 nM of primers, 2.5 [A of cDNA (1:10 
dilution) and the total volume made up with sterile- 
distilled water. The cycling conditions for all primer sets 
consisted of 95°C for 15 minutes followed by 40 cycles 
of 60°C for 10s, 95°C for 30s, 95°C for 1 minute, 55°C 
for 30s and 95°C for 30s. Three biological replicates were 
run for each sample on a plate. Each gene was measured 
in at least six real-time PCR runs (replicate reactions/ 
biological repeats). The real-time PCR (qPCR) reactions 
were run along with two normalising control' genes, 
ribosomal S7 and elongation factor EF1. Primer pairs 
exhibiting a linear relationship between Ct values and 
template concentration in standard curves, and high 
PCR amplification efficiency were chosen for further 
analysis. 

Data analysis of gene expression 

Outlier data points were identified and excluded from 
analysis based on obvious deviations in both the normal 
shape of amplification curves and the Ct values of other 
repeated observations (biological triplicates). The mean 
of the threshold cycles (Ct) for each gene were normal- 
ised against the average values for ribosomal S7 and 
elongation factor EF1. To compare gene expression be- 
tween treatments, we calculated the AACt [51]. The 
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fold-change in gene expression for each target gene, 
normalized to the ribosomal S7 and elongation factor 
EF1 relative to a susceptible strain of An. arabiensis 
from Dar es Salaam, was calculated according to the 2" AACt 
method incorporating PCR efficiency [51]. Basic data ana- 
lysis (regression and t-tests were performed in Excel, with 
p < 0.05 used to assess significant difference between treat- 
ments for the t-tests. 

Ethical considerations 

The study received ethical clearance from the KCMC 
Research Ethics Committee. 

Results 

Microarray experiments involving Moshi samples indi- 
cated that 87 probes were significantly up regulated in 
all three comparisons and 85 down- regulated coming 
from a total of 133 genes (Figure 2). There was a high 
correlation among fold changes between the different 
experimental comparisons, suggesting consistency of re- 
sults (Figure 3). Surprisingly, the genes that were down 
regulated include multiple detoxification genes, includ- 
ing several that have come up as top upregulated hits 
in previous An. gambiae s.s. arrays, including P450s re- 
peatedly resistance-associated in previous work on An. 
gambiae s.s. (CYP6M2 and CYP6P3). Among the up- 
regulated genes were 3 independent ABC transporters, 
one of which ABCB4 (represented by 3 alternate tran- 
scripts) exhibited almost all significant probes. Also 3 
probes for GSTe7 were up-regulated. Full microarray 
results from both experiments have been submitted as 
Additional file 1. However, few genes were selected for 
qPCR. Of about thirty primers that were designed, only 
six passed the criteria for qPCR. These primers are 
shown in Table 1 below including the primers for 



CYP4G16, the gene that was documented in other 
studies reporting on insecticide resistance in An. ara- 
biensis and two normalising control' genes, ribosomal 
S7 and elongation factor EF1 (Table 1). Only two genes, 
ABC 2060 and CYP4G16 were the most significantly 
over-transcribed genes in qPCR (Figure 4). 

Discussion 

Transcriptional analysis of pyrethroid resistance in the 
four day old Lower Moshi An. arabiensis females was in- 
vestigated using the An. gambiae whole genome micro- 
array [25] and quantitative real-time qPCR. Since a fully 
sequenced genome for An. arabiensis was not available; 
we used heterologous hybridization of the An. arabiensis 
samples to an An. gambiae microarray. The microarray 
analysis of this study clearly indicates that the differ- 
ences among arrays could not be appropriately deter- 
mined. In the previous studies, the use of An. gambiae 
'detox chip' could successfully determine the gene ex- 
pression of pyrethroid resistant populations of other spe- 
cies including An. funestus [53], An. stephensi [54] and 
An. arabiensis [28]. Gene expression of pyrethroid resist- 
ant An. arabiensis populations from Pemba Island was 
also successfully determined by using An. gambiae whole 
genome microarray [40]. The microarray analysis of this 
study clearly indicates that there were some very inter- 
esting hits among the up regulated genes most notably 
three independent ABC transporters and GSTe7. 

Transcription of four genes was evaluated in this study 
using quantitative real-time (qPCR). It included one ABC 
transporter gene ABC 2060 that was up regulated, two 
P450s (CYP6M2, CYP6P3) that were down regulated in 
microarray experiments but previously associated with 
pyrethroid resistance and one P450 gene CYP4G16 that 
was identified through published literature and implicated 
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Figure 2 Volcano plot of all probes significant in all of three experimental comparisons of An. arabiensis permethrin-exposed survivors 
vs. unexposed samples (average of the three experiments). Red triangles show probes of candidate genes chosen for subsequent qRT-PCR, 
with probes for the same gene within dashed ovals; the tight cluster of overexpressed gene probes are grouped together. 
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in permethrin resistance in An. arabiensis [40]. Only two 
genes, CYP4G16 and ABC 2060, showed high expression 
levels (p < 0.05) in the pyrethroid resistant Lower Moshi 
An. arabiensis and it is possible that the two genes may be 
playing a role in the observed resistance to pyrethroids 
in the Lower Moshi An. arabiensis population. The low 
levels of expression of the two P450s, CYP6M2 and 
CYP6P3 in the pyrethroid resistant Lower Moshi An. 
arabiensis implies that the observed pyrethroid resist- 
ance results from other mechanisms rather than in- 
secticide metabolism. 

CYP4G16 has previously been recorded as an over- 
transcribed gene in An. arabiensis strain from Cameroon 
[28], Pemba Island Zanzibar [40] and in a laboratory strain 
of An. arabiensis from Sudan that was resistant to DDT, 
permethrin and deltamethrin [41]. In Pemba Island up- 



regulation of CYP4G16 in resistant An. arabiensis ranged 
from 2.0 to 4.5 fold compared to susceptible strains of 
An. arabiensis from Dar es Salaam and Unguja and it 
has been suggested that it plays a role in cuticular-based re- 
sistance. Cytochrome P450s CYP4G17 and CYP4G16 of 
the Anopheles gambiae are closely related to the Drosophila 
melanogaster CYP4G1, the specific insect P450 with decar- 
bonylase function that is involved in hydrocarbon biosyn- 
thesis. CYP4G16 and CYP4G17 proteins are found in the 
mosquito head and the carcass of the abdomen, MozAtlas 
[55] but not the midgut or malpighiam tubules, where 
P450 detoxification genes have been localized previously. 

In this study, metabolic resistance to pyrethroids in 
Lower Moshi An. arabiensis has been shown to be associ- 
ated with increased transcription of P450 gene, CYP4G16 
and ABC transporter. Following qPCR analysis, CYP4G16 



Table 1 Primers used in quantitative real-time PCR (qPCR) (F 


= forward, R = reverse) 




Primer 


Primer sequence 


Fragment length (bp) 


Reference 


ABC 2060 


F: 5'-AATGCACTGCmGCGAACT- 3' 
R: 5'-GACCATCCCACTGmCGGA- 3' 


90 


This paper 


CYP6P3 


F: 5 '-GTG A^G ACG A A ACCC^CGG A AGT-3 ' 
R: 5'^CACCAGTGTTCG<TTCGGGA-3' 


N/A 


Witzig et al, 2013 [52] 


CYP6M2 


F: 5'-TACGATGACAACAAGGGCAAG- 3' 
R: 5'- GCGATCGTGGAAGTACTGG-3' 


N/A 


Witzig et al, 2013 [52] 


CYP4G16 


F: 5'- AGCTG A ACGG ATACCTGG ACCG A-3 ' 
R: 5'- AACACGGAGTGCGAACTGCCAAC-3' 


81 


Jones et al, 2013 [40] 


S7 


F: 5'- AGAACCAGCAGACCACCATC-3' 
R: 5'- GCTGCAAAOTCGGCTATO-3' 


149 


Jones et al, 2013 [40] 


EF 


F: 5'- GGCAAGAGGCATAACGATCAATGCG-3' 
R: 5'- GTCCATCTGCGACGCTCCGG-3' 


N/A 


Jones et al, 2013 [40] 
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Figure 4 Relative expression levels of four candidate genes between sample groups of An. arabiensis differing in predicted susceptibility. 

Data are delta-delta CT values calculated from 5-9 biological repeats. **P < 0.01; *P < 0.05 (t-tests). 



was over-expressed in Lower Moshi pyrethroid resistant 
An. arabiensis compared to susceptible strains of An. ara- 
biensis from Dar es Salaam. (FC = 3.92, p = 0.007). The 
gene ABC 2060 was over-expressed with the FC of 5.07, 
(p = 0.009). ABC transporters have been linked to insecti- 
cide resistance in several species such as Helicoverpa 
armigera [56], Aedes spp. [57,58], Culex pipiens complex 
[59] and Anopheles gambiae [60]. ABC transporters have 
also been correlated with insect resistance to Bacillus 
thuringiensis [61]. Members of the ABC-transporter 
family (ABC-transporters) pump foreign molecules out 
of insect cells using an ATP-dependent mechanism 
[62-65]. They are also involved in lipid transport to the 
cuticle [62]. Interestingly, AGAP002060 is a member of 
the ABCH family, which have been implicated in trans- 
portation of lipids to the cuticle [62]. The previous 
study on An.arabiensis populations a few kilometers 
from Mabogini village [45] (the current study site) re- 
vealed over-transcription of esterase AGAP006227, the 
UDPGT AGAP006775 and the cytochrome P450s CYP9J4 
and CYP6P1. However, its involvement in insecticide re- 
sistance is not clear. In addition, the gene AGAP000987 
encoding the cuticle protein CPAP3-Alb was also 
strongly over- transcribed in An. arabiensis in the area. 
Cuticle proteins play an important role in pyrethroid 
resistance through altered insecticide penetration and 
over-transcription of genes encoding them have been re- 
ported in pyrethroid resistant species including Anopheles 
species [45,66]. 

Conclusions 

The study compared the gene expression profiles of 
permethrin-resistant and susceptible populations of An. 



arabiensis using quantitative real-time PCR (qPCR). 
CYP6P3 is one of the leading candidate genes for meta- 
bolic resistance in West African An. gambiae and its down 
regulation suggests that quite different mechanisms may 
operate in the different species. 

Based on the quantitative real-time PCR q(PCR), one 
P450 monooxygenase gene CYP4G16 and ABC trans- 
porter ABC 2060 were over-transcribed in permethrin 
resistance in Lower Moshi An. arabiensis. However, fur- 
ther investigation is needed to clarify its role in observed 
resistance. There is no real evidence that ABC 2060 
transporter gene is linked to observed resistance, al- 
though significant lower expression of this gene was ob- 
served in a susceptible strain from Dar es Salaam than 
the Moshi population (i.e. either survivors or Moshi 
control). However, the two genes identified are most 
likely not the only genes involved in pyrethroid resistant 
An. arabiensis from Lower Moshi. Further investigation 
is needed especially in microarray analysis to identify 
more genes that may be involved in insecticide resist- 
ance in the country and its effects on the current mal- 
aria control interventions. 

Additional file 



Additional file 1: Full microarray results from both experiments. 
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